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We study corrections to the conformal hyperscaling relation in the conformal window of the large 
Nf QCD by using the ladder Schwinger-Dyson (SD) equation as a concrete dynamical model. From 
the analytical expression of the solution of the ladder SD equation, we identify the form of the leading 
mass correction to the hyperscaling relation. We find that the anomalous dimension, when identified 
through the hyperscaling relation neglecting these corrections, yields a value substantially lower than 
the one at the fixed point 7^ for large mass region. We further study finite- volume effects on the 
hyperscaling relation, based on the ladder SD equation in a finite space-time with the periodic 
boundary condition. We find that the finite-volume corrections on the hyperscaling relation are 
negligible compared with the mass correction. The anomalous dimension, when identified through 
the finite-size hyperscaling relation neglecting the mass corrections as is often done in the lattice 
analyses, yields almost the same value as that in the case of the infinite space-time neglecting the 
mass correction, i.e., a substantially lower value than 7* t for large mass. We also apply the finite- 
volume SD equation to the chiral-symmetry-breaking phase and find that when the theory is close 
to the critical point such that the dynamically generated mass is much smaller than the explicit 
breaking mass, the finite-size hyperscaling relation is still operative. We also suggest a concrete 
form of the modification of the finite-size hyperscaling relation by including the mass correction, 
which may be useful to analyze the lattice data. 
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I. INTRODUCTION 



Technicolor model 0, 0] has been considered as an interesting possibility for the dynamical origin of the electroweak 
symmetry breaking. However, it has fatal phenomenological difficulties (especially with the strong suppression of 
{y-^ I flavor changing neutral current processes). The problems can be solved by the walking technicolor 0, Q having 
approximate scale invariance with large mass anomalous dimension, 7 m ~ 1 , which was proposed based on the ladder 
Schwinger-Dyson (SD) equation. Modern technicolor models often utilize asymptotically free gauge theories with an 
approximate infrared fixed point (IRFP) to achieve the walking behavior. 

The SU(iV) gauge theory with a large number of massless fermions is one of the theories that are expected to possess 
£NJ \ such a property [5]. In the case of SU(3) gauge theory for example, the two- loop running coupling has an IRFP in 
the range of 9 < Nf < 16 (< N^), where Nf is the number of massless fermion with fundamental representation, 
and N^ F is the value of Nf above which a theory loses its asymptotic freedom nature @, 0] ■ Within this range of 
Nf, the larger the number of Nf becomes, the smaller does the value of the running coupling at the IRFP. Because 
of this, it is expected that there is a critical value of flavor, N", below which the theory is in the confining hadronic 
phase with broken chiral symmetry, while above which it is in the deconfined phase with unbroken chiral symmetry. 
An analysis based on the SD equation with the improved ladder approximation estimates that the value of Nj 1 lies 
between 11 and 12 Therefore, for 12 < Nf < 16 (often called "conformal window"), the theory possesses an exact 
IRFP, while for 9 < Nf < 11, the chiral symmetry is spontaneously broken, i.e., the IRFP disappears and the scale 
invariance is only approximate. In Ref. [a. H(ij|. this chiral phase transition at N" was further identified with the 
"conformal phase transition" which was characterized by the essential singularity scaling (Miransky scaling) . 

Considering the intrinsically non-perturbative nature of the problem, the lattice gauge theory should play an 
important role for the study of the phase structure of such theories. In addition to pioneering works such as Refs. [TTI - 
[3], there is growing interest in this subject in recent years [To] ]. A straightforward way of investigating the infrared 
behavior of a given theory is to calculate the running coupling constant of the theory. Though it requires simulations 
in a wide range of parameter space since extensive range of the energy scale has to be covered to trace the running of 
the coupling by step-scaling procedure, there are many groups that devote their efforts to such a direction. 

Alternatively, infrared conformality of the theory can also be investigated by deforming the theory with the in- 
troduction of a small fermion bare mass, too, as a probe, and study the relation between some low-energy physical 
quantities (such as the meson masses and the decay constants) and mo. In Ref. |jil[lj}, it is shown that the scaling 
relation between a low-energy quantity and toq can be expressed in terms of the mass anomalous dimension at the 
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IRFP, 7* n . 1 In the case of the mass (M) of a meson with certain spin and quantum numbers, for example, the scaling 
relation ( "hyperscaling relation" ) is expressed as 

M~mJ /(1+ ^. (1) 

When one considers a theory in a finite space-time, the scaling relation is modified to the "finite-size hyperscaling 
relation" as follows: 

M = L- 1 f(x), (2) 

where, L is the size of space and time, and / is some function of scaling variable x which is defined as 

x = LmJ /(1+7 ™ ) . (3) 

Here, we introduced dimensionless quantities, L = LA and mo = mo/ A, where we take A as the UV scale at which the 
infrared conformality terminates. Several groups [18l - l21| tried to judge whether candidate theories posses an IRFP or 
not by measuring the low-energy quantities on the lattice for various combination of input values of L and mo, then 
checking whether Eq. @ is satisfied for a certain value of ^* n . 

However, a couple of questions arise here regarding use of (finite-size) hyperscaling relation for the study of infrared 
conformality: One of them is related to the fact that the bare fermion mass, mo, which is introduced as a probe, itself 
necessarily breaks the infrared conformality of the original theory. How small mo has to be so that the hyperscaling 
relation is approximately satisfied? What is the form of correction if it is not small enough? When the anomalous 
dimension is measured for mass not so small, can it be regarded as 7* n at IR fixed point at face value? Another 
question is, when the theory in question does not have an IRFP (namely, in the phase where the chiral symmetry is 
spontaneously broken) in the first place, how and how much is the hyperscaling relation violated? 

The ladder SD equation, which is the birth place of the walking technicolor, is actually a concrete dynamical model 
to study such questions. In the framework of the ladder SD equation, we know whether a given theory is infrared 
conformal or not, and also the value of the anomalous dimension as well. Therefore, we can quantitatively study 
the (finite-size) hyperscaling relation and its violation by using the ladder SD equation in a self-consistent manner. 
Numerical calculations can be easily done in a wide range of parameter space, and to a certain extent, even an 
analytical understanding can be obtained by investigating the solution of the ladder SD equation. 

In this paper, we study the (finite-size) hyperscaling relation and its violation, based on the ladder SD equation by 
taking the example of SU(3) gauge theory with various number of fundamental fermions (which is often called the 
large N f QCD). 

In the next section, from the analytical expression of the solution of the ladder SD equation, we identify the form 
of the leading correction to the hyperscaling relation. We find that the anomalous dimension 7 m for off the IRFP 
(with finite mass scale) is substantially smaller than 7*„, the value at IR fixed point (with vanishing mass) for the 
larger mass. Our result may shed some light on the value of the anomalous dimension often reported by the lattice 
simulations done with relatively large masses. 

In Section 3, for the purpose of studying the finite-size hyperscaling relation, we formulate the SD equation in 
a finite space-time with the periodic boundary condition. By numerically solving it for various values of the input 
parameters (jho,L), finite-size, as well as mass deformation effects on the finite-size hyperscaling relation is studied 
in the conformal window. The result suggests that the correction due to the finite-size effect on the finite-size 
hyperscaling relation is negligible compared with that coming from large mass corrections. Then we find that the 
anomalous dimension, when identified through the finite-size hyperscaling relation neglecting the mass corrections, is 
almost the same as that obtained through the hyperscaling in the infinite space-time neglecting the mass corrections 
and hence is substantially smaller than 7^. We also use the ladder SD equation in a finite space-time to study the 
chiral-symmetry-breaking phase, and show how the effect of spontaneous breaking of the chiral symmetry affects the 
(finite-size) hyperscaling relation. In the case of Nf = 11 which is close to the criticality so that the dynamically 
generated mass is small compared with the explicit mass, the finite-size hyperscaling relation is still operative. We 
further suggest a concrete form of the modification of the finite-size hyperscaling relation taking account of the mass 
correction, which may be useful for the analysis of the lattice data. 

Finally, Section 5 concludes the paper. 



This 7* n is identified with the anomalous dimension f m relevant to the walking technicolor, which is ■y T n measured at ultraviolet (UV) 
limit (instead of IR limit), or near the scale of (pseudo-) UV fixed point, usually identified with the ETC scale. See discussions below 
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FIG. 1: Two-loop running coupling (solid curve) compared with the approximate form in Eq. ([$]) (dashed line) in the case of 
SU(3) gauge theory with 12 massless fundamental fermions. 





12 


13 


14 


15 


16 


a* 


0.75 


0.47 


0.28 


0.14 


0.042 


7m 


0.80 


0.36 


0.20 


0.095 


0.027 



TABLE I: Values of a, and also listed are the corresponding 7^ for SU(3) gauge theory with fundamental fermions in the 
conformal window to be given by Eq. (|28[) . 



II. CORRECTION TO THE HYPERSCALING RELATION 



We start from the study of the hyperscaling relation in the infinite space-time, namely the one in Eq. (p}. In this 
section, from the analytical expression of the solution of the ladder SD equation, we identify the form of the leading 
correction to the hyperscaling relation. Here, we take the example of the large Nf QCD in the conformal window to 
study infrared conformal theories. 

The two-loop running coupling of the large Nf QCD is shown in Fig.[TJ In the figure, the two-loop running coupling 
(and its approximated form) in the case of SU(3) gauge theory with 12 massless fundamental fermion is plotted as an 
example. The solid curve represents the two-loop running coupling which are obtained from the following RGE for 

«(M 2 )(=^¥ 



a( M 2 ) = p{a{^)) = -ba 2 ^ 2 ) - ca 3 (» 2 ), 
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a* is the value of the running coupling at the IRFP which is determined as 

b 



N c 
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(4) 



(5) 



(6) 



Values of a* in the case of SU(3) gauge theories with various number of fundamental fermion (in the conformal window) 
are shown in Table |TJ A which appears in Fig. Q] is a renormalization group invariant scale, A = \x exp ^— J a ^ -^^j 1 
the two-loop analogue of Aqcd of the ordinary QCD, which is taken as Q 



A = \i exp 



b a. .. 



■log 
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b a(fj?) 



a(A 2 ) ~ 0.78a*, 



(7) 
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in such a way that the scale A plays the role of the "UV cutoff " where the infrared conformality we are interested 
in terminates, i.e., a(/i 2 ) ~ const(~ a„) for /i 2 < A 2 , while a(p, 2 ) ~ 1/ log(/i 2 /A 2 ) for /.t 2 > A 2 as in the usual 
asymptotically free theory. Actually in the walking technicolor, A is taken to be of order of (or even larger than) 
the UV scale Aetc ("ETC scale") where the technicolor theory no longer makes sense as it stands and is actually 
converted into a more fundamental theory such as the Extended Technicolor (ETC). 

In this paper we use the following form of the running coupling as an approximation of the two-loop running 
coupling (dashed line in Fig. [1]) : 

In this approximation the coupling takes the constant value a* (the value at the IR fixed point) below the scale A 
and entirely vanishes in the energy region above this scale. Therefore, the physical picture of the large Nf QCD with 
this approximation is the same as that in constant coupling gauge theory with UV cutoff A which was extensively 
studied long time ago [22| - |24j . We note that it is possible, at least numerically, to solve the SD equation without this 
approximation for the two-loop running coupling. However, we adopt this simplification so that we can analytically 
study the solution of the SD equation to a certain extent. 



A. SD equation 

Let us first write down the SD equation for SU (N c ) gauge theory with fundamental fermions: 

iS-/(p) = t - m + J J^L C 2 g 2 ((p - k) 2 ) ^-L_ ^ _ {p -^} P k - 2 k) » ) >f iS F {k) 7 *. (9) 

Here, iSp 1 = A(p 2 )p' — B{p 2 ) is the full fermion propagator, C% — is the quadratic Casimir, and g((p — q) 2 ) is 

the running coupling constant. In the above expression, we took the Landau gauge, and adopted the improved ladder 
approximation, in which the full gauge boson propagator is replaced by the bare one and the full vertex function is 
replaced by a simple 7^-type vertex with the running coupling constant associated with it. From this equation, we 
obtain the following two independent equations: 



i( 2, _ , , f d 4 k C 2 g 2 ((p-k) 2 ) 2 [ (P' k ) , Jp-(p-k)}{k-(p-k)} 

A(P ) — 1+ ——J TTT7T2T2 i d7l2Y2 ^ > 



(2tt) 4 k 2 A(k 2 ) 2 + B(k 2 ) 2 ' 



2 

p 2 (p — k) 2 p 2 (p—k) A 



2 . _ , d 4 fc 3C 2 g 2 ((p-k) 2 ) B(k 2 ) 



(10) 

B(p z ) = m + J -j—^ fc2j 4 ( p)2 + B(fc 2 ) 2 (P - k) 2 ' (U) 

These are coupled equations for A(p 2 ) and B(p 2 ) written in Euclidean momentum space. (Note that we have dropped 
the subscript "E" for Euclidean momentum variables.) 

To further simplify the SD equation, we adopt a simplified form for the argument of the running coupling, g ((p — 
k) 2 ): we take it to be a function of only p 2 and k 2 instead of (p — k) 2 . With this simplification, it becomes possible to 
carry out the angular integration (in the momentum space), then the SD equation becomes an equation for a single 
variable x = p 2 E . Also, in this case, A(x) = 1 is obtained from Eq. (|10l) . Therefore, the SD equation becomes a single 
integral equation for the mass function £(#)(= B(x)/A(x) = B(x)). 

For the analytical study we adopt a practically simple ansatz for the running coupling: pH ] 

-g 2 ({p-k) 2 )^g 2 {max{p 2 ,k 2 }). (12) 
Then the ladder SD equation with Eq.® reads: 

w ^ x 3 ° 2 t\ 1 ^ 

T,{x) = rn / dy , (13) 

4?r J q max{x,y\ y + Y,^(y) 

which can readily be converted into an equivalent (nonlinear) differential equation with boundary conditions (23|: 

+ <"§OT " °' (14) 

lima; 2 E(x)' = 0, (15) 
(xS(x))'\ . a = too . (16) 
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We may further simplify Eq. (fl4l) by replacing S(x) in the denominator of the second term in the LHS by a constant, 
mp, which is customarily defined by : 

rap = E(.t = m P ). (17) 

Then the SD equation reads: [26[ 

( x Z {x) f + am ^JM^ = o. (18) 
An x + mp 

We should note that it is known that the solution obtained from this linearized equation well approximates that 
obtained (by numerical calculation) from the equation without linearization. Later we shall show that particularly for 
the anomalous dimension there is a remarkable agreement between the analytical result obtained from the asymptotic 
solution of the linearized equation Eg. (1181) and the numerical one from the full nonlinear integral SD equation under 
a slightly different ansatz for the argument of the running coupling (to be mentioned later). 

m p defined in Eq. (JT7J) is often called the "pole mass" , though of course mp is not the real pole mass (remember 
that x is the Euclidean momentum square). Still, mp is useful quantity for the investigation of the hyperscaling 
relation since it is known, from the study with the Bethe-Salpeter equation (27| . that mp is proportional to meson 
masses. Therefore, in this paper, we use mp which is obtained from the solution of the SD equation as low-energy 
physical quantity which appears in the hyperscaling relation. 



B. Leading correction to the hyperscaling relation 

Before we proceed to the investigation of the solution of the SD equation to derive the relation between mp and too, 
we note that it is known [22J that there is a critical value of a*, a cr (^ 0), such that spontaneous symmetry breaking 
solution ^ 0) which satisfies Eqs. (|18p. (TT51) and p^|) for the chiral limit too = mo(A) = does not exists for 



where [2S 



(19) 



acl = oTT = T ( N c = 3) ■ (20) 



3C; 



2 

Namely, a nontrivial solution (S(a;) ^ 0) for a < a CI is the explicit breaking solution which exists only for to — 
too (A) 7^ 0. Then the pole mass mp is nothing but a renormalized mass ("current mass") mp: 

m P = m R = Z^ to , (21) 

where mp = mp(p = mp), and Z m = Z m \^ =mp is the mass renormalization constant. This means that, in the 
chiral limit mp = 0, there is no mass gap (dynamically generated mass) mp, — 0, and therefore the IRFP of the 
theory is exact for a* < a cr . This is exactly the region in which one expects that the hyperscaling relation should be 
satisfied. Therefore, in the rest of this section, we concentrate on studying the SD equation in the region of a* < a cr , 
or equivalently (through the relation in Eq. ©) > Nf > N? (conformal window), where = 16.5 and 

iV^ r ~ 11.9 for N c = 3. 

A solution of Eq. (|18l) which satisfies boundary condition Eq. (|15p can be expressed in terms of the hypergeometric 
function as [26] 

Z( X )=(rn pF (^^,2*), (22) 



"p 



where 



W = J1-^L. (23) 



£ is a numerical coefficient which is determined from the definition of mp in Eq. (1171) : 



1^,1^,2,-1). (24) 
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In the limit of j:> mp, the solution can be expanded as 



Ti{x) ~ £mp 



r(«) 



r(^)rW) 



— J" + (w O -w) 



By inserting the above expression of into the remaining boundary condition in Eq. (H 
relation between mp and mo: 



(25) 



, we obtain the following 



m 



£m P 



r(H±l)a 



A 2 



+ (wf> —a;) 



From this we have the mass renormalization constant Z m in Eq. (|21|) as 



(26) 



m _ mo_ _ 
% m P 



r( w ) 



T(^)2 ^ mp 



A 2 . 

IT + (W ^ -W) 



(27) 



where we note again mp = mp in the conformal window. 

Then we can obtain the mass anomalous dimension at IRFP, 7* 



as 



7„ 



lim 



5 log Z n 



m P /A->o 91og(mp/A) 



= 1-u = I-4/I- 



(28) 



where the limit is taken as mp — >• with A fixed. In Table U we show values of 7^, which can be calculated from 
the above expression combined with Eqs. (|5j), JBJl and (l20l) . in the case of SU(3) gauge theory with various numbers 
of fundamental fcrmion in the conformal window. 

It should be noted that this 7^ at IRFP is actually the same as the anomalous dimension in the UV limit (A — > 00 
with mp fixed) 7 



(UV) 



lim 



with tJ? = 1 (the value at (pseudo-) UV fixed point in the broken phase a* > a CI in the chiral limit mp = 0: mp = 
m D ) 3J: The technifermion condensate (^V)U at the UV scale A = A ET c(> 10 3 TeV) » fi (= 0{m P ) = O(TcV)) is 

enhanced by 7 i UV) = 1 as (^)U = ^ m 1 (^)U with z Vn = (A//*) 7 ™^ = (A/m) 1 » 1. 
By using this 7^, we can rewrite the expression in Eq. (|2"B"|) in terms of 7^: 



d log Z m 
A/mp^oo aiog(mp/A) 



1 — uj 28] . 7m is the quantity relevant to the walking technicolor 



m _ 
A 



r(2^>)2 V a / 



r(^) 2 



V A 



3-7; 



(29) 



This is the expression which should be compared with the hyperscaling relation in Eq. ([T]). It is obvious that if we 
drop the second term in the RHS of Eq. (|2Q[) . it reduces to the hyperscaling relation [lq |. Therefore, the second term 
should be identified as the leading correction to the hyperscaling relation. 

To see the significance of the correction term, in Fig. [2J we plot ratios of the second term to the first term in the 
RHS of Eq. ([29]) as functions of mp/A for various values of 7^ in the range of < 7^ < 1.0, which corresponds to 
N^ F > Nf > Nf X . When mp is much smaller than A, (except in the case of 7* t = 1.0) the effect of the second term 
is very small since in the range of < 7*„ < 1.0, the power of (mp/A) in the second term is always greater than 
that in the first term. This is reasonable considering the fact that small mp (or equivalently, small mo) means small 
mass deformation. Another limit in which Eq. (|29p approximates well the hyperscaling is 7*„ — > 0. In this limit, the 
coefficient of the second term goes to while that of the first term goes to 1. Also, the power suppression of the second 
term becomes strong in this limit as well. However, we should remember that phenomenologically motivated theories 
have large 7^. In the limit of 7^ — > 1, the power of (mp/A), as well as coefficients of the two terms asymptote to 
the same values. Therefore, we have to take the second term seriously when we study the anomalous dimension of 
the candidate theories for viable walking technicolor models with 7 m ~ 1 through the hyperscaling relation from the 
numerical data on the lattice. 

For checking the reliability of our ansatz Eq. (fT2l and the linearization of the differential equation, as well as the 
asymptotic expansion Eq. (|25"]) . we show in Fig. [3] the log-scale plot of mo — mp in Eq. ([29]) for Nf = 12 in comparison 
with that obtained by directly solving numerically the full nonlinear SD equation with more natural (angle-averaged) 
ansatz: 



g 2 ((p-fc) 2 )^gV + fc 2 ) 



(30) 
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FIG. 2: Ratios of the second term to the first term in the RHS of Eq. (|29[) as functions of mp/A for various values of 7^ 
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Mp/A 



FIG. 3: Analytical asymptotic solution of the linearized SD equation vs. numerical solution of the full nonlinear one (for 
Nf = 12). Analytic result (blue dotted line) is the plot of Eq. ()29[) which is from the asymptotic expansion of the linearized 
SD equation Eg. (|18[) with the ansatz Eg. (|12[l . while the numerical one (denoted by red plus) is that of the solution of the full 
nonlinear integral SD equation Eq. (|31[) with ansatz Eq. (|30[) , 



The SD equation in this case with Eq.© reads: 

S(x) = too + / dy , (31) 

4?r J max{x,y\ y + r, z (y) 

which differs from Eq.(JT3J) by A 2 — x in the UV end of the integral. The slope in Fig.© corresponds to 7,„ + 1. The 
agreement is remarkable, irrespectively of the different ansatz and the additional approximations. 
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FIG. 4: Effective mass anomalous dimension as a function of mp/A for SU(3) gauge theories with 12, 13, 14, 15 and 16 
fundamental fermions. 



C. Effective anomalous dimension 



In the literature, the hyperscaling relation is often used as a tool to judge whether a theory is infrared conformal 
or not. However, the importance of the corrections to the hyperscaling relation due to the mass deformation are 
often underestimated, or even completely neglected. This is not surprising because, in practical situations, it is very 
difficult to notice that data (for example, hadron mass, Mh, obtained from lattice simulations with various values of 
input moa) need to be fitted by a function with correction term. Even in the situation that the correction term is not 
very small compared to the leading term, data could be easily fitted by a simple function of a form m ~ unless 
data are taken in a wide range of mo (or, equivalently, My). Especially when the data are associated with, say, a few 
percent of error bars, it is very possible that one succeeds in fitting the data with a function of a form uiq ~ M^ +7 . 
However, the best-fit value of 7 obtained by this fitting must be numerically different from the actual mass anomalous 
dimension at the IRFP. To make the difference clear, we introduce effective mass anomalous dimension, 7^, which is 
defined as the value of 7 one obtains as a best-fit value when one forces to do fitting by using a fit- function which has 
a form of hyperscaling relation. Since the significance of the correction term is different for different values of Mb, 
the value of 7^f should change depending on the range of Mh one uses for fitting to obtain it. 

In the framework of the SD equation with the improved ladder approximation, we can identify 7^f as 7 m = 
7 m (/^/A)| M=mp obtained from Eq.([2"Tf with Ea.([28"]l (or equivalently Eq. ((29)) ): 

ff d\ogZ m 

7m —7m — fl1 „„/ /A \ l-^J 



91og(m P /A) 
d 

aiog(m P /A) 



log k 



r(l-7») / m P y"m r(-l + 7*) ( m P y-ri 
r(2z_>)2 V a / r(^) 2 v a J 



(33) 



It is obvious that were it not for the second term, 7^ would coincide with 7 m . Note again that the significance of the 
correction term is different for different values of mp. Therefore, the effective mass anomalous dimension becomes a 
function of mp. In Fig.|H 7^ for SU(3) gauge theories with 12, 13, 14, 15 and 16 fundamental fermions are plotted as 
a function of mp. For the purpose of making it easier to see the deviation of the effective mass anomalous dimension 
from the value at the IRFP, we also plot 7^/7^ in Fig. [3] From these figures, as we expected, we see that the 
deviation between 7^ and 7 m becomes more significant in larger mp/A region. We can also see that the deviation of 
the effective mass anomalous dimension is larger for smaller Nf (or, in other words, for Nf closer to N" ~ 11.9). This 
is also expected from the discussion below Eq. (|2T)1) because smaller Nf means larger 7* n , with which the correction 
term to the hyperscaling relation becomes important. As we mentioned earlier, phenomenologically interesting theory 
is the one with large mass anomalous dimension. Therefore, it is important to keep this effect of correction term to 
the hyperscaling relation in mind when one study such theories. 
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FIG. 5: Effective mass anomalous dimension which is normalized by the value of it at the IRFP as a function of mp/A for 
SU(3) gauge theories with 12, 13, 14, 15 and 16 fundamental fermions. 



III. EFFECT OF THE CORRECTIONS ON THE FINITE-SIZE HYPERSCALING 



So far we have concentrated on the study of the hyperscaling relation in the infinite space-time. However, since the 
lattice simulations are done in a finite space-time, the hyperscaling relation in a form of Eq. ([2]) are used more often. 
Therefore, it is important to study the effect of correction due to mass deformation on the finite-size hyperscaling 
relation. For the purpose of studying the finite-size hyperscaling relation, we formulate the SD equation in a finite 
space-time with the periodic boundary condition. By numerically solving it for various values of input parameters 
(mo,I), mass deformation effects on the finite-size hyperscaling relation is studied in large Nf QCD. 



A. SD equation in a finite space-time 



To formulate the SD equation in a finite space-time, we start from the SD equation in the infinite space-time in 
Eqs. (fTUf and (fTTjl . To put these equations in a finite space-time, all one needs to do is to replace the continuum 
momentum variables by the discrete ones: 

Pi^Pi = ^ 1 , KeN), (34) 

where, pi is the i-th component of the momentum variable. We adopted the periodic boundary condition for all 
directions, though it is easy to implement the anti-periodic boundary condition. We also assumed that the size of 
all space-time directions are the same. It is also straightforward to introduce different sizes for spacial and temporal 
directions. However, we took the same length for every direction just for simplicity, n^s are integers which label 
discrete momentum variables. With this replacement of the momentum variables, the SD equation in Eqs. (|10p and 
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(fTTj) turn into the following form 2 

1 



A(p) = 1 
B(p) 



^ k 2 A(kf 

"10,1,2,3 V ' 



B(k) 



-A(k) 



(p-k) {p ■ (P - ~k)} {k ■ (P - ~k)} 



p 2 (p — k) 2 



p 2 (p-k) A 



m 



1 

IA 



E 

™0,1,2,3 



3C 2 g 2 ((p-k) 2 ) B(k) 



k 2 A(k) 2 + B(k) 2 (p-k) 2 



(35) 
(36) 



where 



2tt 

T 



/ n \ 

n 2 

w 



2tt 

T 



/ to \ 

TOl 
TO 2 

\m 3 J 



(37) 



Here, n.; and to^ are integers, though we should note that the SD equation is not defined at p — since one of the 
two independent equations is derived by requiring coefficients of j) in LHS and RHS are the same in Eq. (J9j) . In the 
above expressions, A(p 2 ) and B(p 2 ) were replaced by A(p) and B(p). This is because they are no longer functions of 
momentum-squared since the rotational symmetry is broken (except some residual discrete rotational symmetry) due 
to the hyper-cubic shape of the finite-size space-time. However, at the practical level, the effect of such rotational- 
symmetry violation is negligible as far as L is taken large enough compared to the scale of relevant physics. This is 
true in the case of the current study. We first numerically solve, by using iteration method, Eqs. (|35p and (131)1) as 
coupled equations for A(p) and B(p). Then, to obtain the value of mp (see Eq. (j!7p ). we plot £(i?) = B(p)/A(p) as 
a function of p 2 . There are always multiple p's which give the same value of p 2 , and those do not necessarily give a 
degenerate value of E unless those are related by the residual discrete rotational symmetry. However, we confirmed 
that such differences are negligible in the region where mp is determined. We should also note that, when we estimate 
the value of mp from Eq. (j!7[) . we used a function which is obtained by interpolating E(p) in momentum space. 
However, we never did extrapolation to the scale below 2ir/L since there is no reliable information below that scale. 
Therefore, we obtain data only when the value of mp is greater (2w/L) 2 . 



B. Finite-size hyperscaling and its corrections 



In this subsection, by numerically solving the finite-volume SD equation formulated in the previous subsection, 
we generate data of mp for various sets of input parameters (LA, mo/ A). We take SU(3) gauge theory with 12 
fundamental fermion as an example here. Then, by using those generated data, we do the analysis based on the 
finite-size hyperscaling in Eq. ([2]). This is a kind of "simulation" of the practical situation we often encounter when 
we study a theory by using data obtained from lattice simulations. An interesting point about doing hyperscaling 
analysis using data generated by the SD equation is that we know that the SU(3) gauge theory with 12 fundamental 
fermions, in the framework of the SD equation, is the infrared conformal theory, and we also know the value of the 
mass anomalous dimension at the IRFP, which is estimated as 7^ ~ 0.80 in this case. (See Table HI) Therefore, we 
clearly see how the finite-size hyperscaling is violated due to the effect of the mass deformation. 

In Fig. [SI we plot the values of mp /A (horizontal axis) for various values of too/A (vertical axis) and LA (indicated 
by different symbols). When we solved the finite volume SD equation, we adopted an angle averaged form (as in 
Eq. pop ) for the argument of the running coupling. As we explained in the previous section, this is the procedure 
which is needed to make the angular integration (in momentum space) possible in the case of infinite-volume SD 
equation, and it is actually not needed for finite-volume SD equation since we numerically solve them by iteration 
without doing angular integration. However, for the purpose of putting finite- and infinite-volume SD equations on 



2 During the summations of rrij, one encounters singularities at rrij = n^. Also, in the limit of B — > 0, the contribution from 
(mo, mi, 7712, m^) = (0,0,0,0) diverges. However, these can be identified as unphysical artifacts considering the fact that these are 
integrablc singularities in the case of infinite space-time. Therefore, in the numerical calculation of the SD equation, we simply drop 
the singular points from summations. (The latter singularity can also be avoided by adopting the anti-periodic boundary condition. We 
did the numerical calculations with the anti-periodic boundary condition in the temporal direction, and compared the solution with the 
one obtained from the periodic boundary condition with the prescription explained above. We found that the difference between two 
are negligible.) 
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FIG. 6: Values of mp/A (horizontal axis) for various values of mo/A (vertical axis) and LA (indicated by different symbols) 
for SU(3) gauge theory with 12 fundamental fermions. Dashed curve is mp as a function of mo which is obtained from the 
numerical solution of the SD equation in the infinite space-time (Eqs. (|10[) and HI}). 



the same ground, we adopted angle averaged argument for finite-volume SD equation as well. We note that when 
we adopt the angle-averaged form for the argument of the running coupling, summations in the finite-volume SD 
equation are restricted in the range of p 2 + k 2 < A 2 . 

In Fig. |6l in each LA, one notices that data are plotted in the range of mp/ A which is larger than a certain value. 
This lower limit comes from the IR cutoff effect which was explained at the end of the previous subsection. Dashed 
curve in the figure is mp as a function of mo which is obtained from the numerical solution of the SD equation in 
the infinite space-time (Eqs. (fTU|) and (fTTT) with the ansatz Eq. (f50|) and Eq.©). In the figure we see that data for 
LA = 20, 25, 30 are almost degenerate, and take values close to the dashed curve for the infinite space-time. This 
means that LA = 20 is large enough that the finite-size effect is negligible for the determination of mp in this mass 
range. 

Now, let us do the finite-size hyperscaling analysis by using data shown in Fig. EI In Fig- El we plot the values oimpL 
as a function of x = LA(mo/A) 1 ^ 1+7 ^ for 7 = 0.2,0.4,0.5,0.6,0.7 and 0.8. If the theory is infrared conformal, and if 
the effect of the mass deformation is negligible, this kind of plot should show good alignment of data when input value 
of 7 is chosen to be 7^, the value of mass anomalous dimension at the IRFP. Here, we know that, in the framework 
of the SD equation, the theory is infrared conformal, and the value of the mass anomalous dimension at the IRFP is 
7^ ~ 0.8. However, plot in Fig. [7] shows no alignment for 7 = 0.8, instead, data are well aligned for 7 = 0.5 and 0.6. 
This suggests that the effect of the correction to the hyperscaling relation due to the large mass deformation appears 
also in the case of finite-size hyperscaling relation. Note that data we used here are in the range of mp/A > 0.3. In 
that range, from Fig. 0] we see that the effective anomalous dimension takes the value 7^ = 0.5 ~ 0.6. This is the 
reason why the data show good finite-size hyperscaling with input value of 7^f = 7 = 0.5 ~ 0.6. We have done the 
same analysis for SU(3) gauge theory with Nf = 14 and 16, and found similar results. 

It is interesting to ask whether there is a function which can be fitted to all the data shown in Fig. [6] We tried the 
following form of fit function, and found that it can be globally fitted to all the data fairly well: 



m (L, m P ) 



C 



Ampil + Bm 2 ^) 1 '^ + - 



1+7 



(38) 



Here, A, B, C and 7 are fit parameters, and it is understood that all the dimensionful quantities are normalized by A. 
This fit function is similar to the form of finite-size hyperscaling relatinon in Eq. @: When one takes B = 0, Lmp is 



expressed by a function of x = Lm, 



1/(1+7) 




The term proportional to B represents the effect of mass correction to the 



hyperscaling relation. The best-fit values we obtained for fit parameters are: A — 1.52, B = — 0.512, C = 0.323 and 
7 = 0.794. It is remarkable that we obtained value of 7 which is quite close to the value of 7^ = 0.8. For comparison, 
we also did fitting with fixing B = 0, and found that the best-fit value of 7 = 0.52. This is consistent with Fig. [71 in 
which it was shown that finite-size scaling (without correction term) is approximately satisfied when 7 = 0.5 ~ 0.6. 
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Y = 0.2 y=0.4 




FIG. 7: Values of mpL obtained from the finite-volume SD equation as a function of x = LA(mo/A) 1 ^ 1+7 - ) for 7 = 
0.2,0.4,0.5,0.6,0.7 and 0.8 in SU(3) gauge theory with 12 fundamental fermions. Data for LA = 12,16,20,25 and 30 arc 
plotted as different symbols. 



Of course, the power of the second term in the RHS of Eq. (|38|) . namely 2 — 27, is specific to the ladder SD analysis, 
though it is worth trying to do fitting lattice data with using the above fit function. One could also make the power of 
the second term in the RHS of Eq. ([551) as free parameter. If the chi-square of the fitting significantly reduces by the 
inclusion of the correction term, or even if the chi-square does not change very much but the value of 7 significantly 
changes, it is very possible that the best-fit value of 7 which is obtained with correction term into consideration is 
close to the actual value of 7^. 
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C. Violation of the hyperscaling relation in theories with spontaneous chiral symmetry breaking 

Here, by the same procedure used in the previous subsection, we study the finite-size hyperscaling relation in theories 
with spontaneous chiral symmetry breaking. In the case of theories with spontaneous chiral symmetry breaking, mass 
gap exists even in the chiral limit, and therefore an IRFP is only approximate. Here, we show two examples: one is 
SU(3) gauge theory with Nf = 9, and the other is that with Nf = 11. The former is an example of a theory which is 
far away from the conformal window, in which the infrared conformality is expected to be largely violated. The latter 
is an example of a theory which resides close to the conformal window, and the breaking of the infrared conformality 
due to the spontaneous chiral symmetry breaking is expected to be small. 

In Fig. IH we show the plots of mpL obtained from the finite-volume SD equation in SU(3) gauge theory with 
9 fundamental fermions as a function of x = iA(mo/A) 1 ^ 1+7 ^ ) for 7 = 0, 0.5, 1.0, 1.5 and 2.0. Data for LA = 
12, 16, 20, 25 and 30 are plotted as different symbols. As we expected, since the infrared conformality is largely broken 
due to the spontaneous chiral symmetry breaking, large violation of hyperscaling relation is observed. Note that the 
dynamically generated mass for Nf — 9 is mp,/ A ~ 0.58, where mp is the value of mp obtained by the spontaneously 
broken solution of the ladder SD equation in the chiral limit too = 0. (mp is estimated by the SD equations, Eq. (|31l) . 
in the chiral limit m = 0.) This is compared with the typical values in Fig. [5J mp/A = 0.58 — 0.77 for LA = 30. 

On the other hand, a similar plot for SU(3) gauge theory with 11 fundamental fermions is given Fig. [§J We show 
the result for 7 = 1.0, with which we found data are best aligned each other. Again, as we expected, since the 
theory is close to the chiral restoration point, and the effect of the spontaneous chiral symmetry breaking is small, 
the violation of hyperscaling relation is small. Note that mp/A ~ 0.05 for Nf — 11, while typical values of mp in 
Fig. [S] are mp/A = 0.28 — 0.69 (3> mp,/ A) for LA — 30. Of course, one can see that there is a small amount of 
misalignment. However, let us imagine those were data obtained from lattice simulations, and each data point has, 
say, a few percent error bar, in which case, the data might look consistent with conformal hyperscaling. Therefore, 
when one obtained data which look consistent with conformal hyperscaling with a large mass anomalous dimension, 
there is a possibility that the theory is exactly the one the technicolor model favors, namely the dynamics with 
spontaneous chiral symmetry breaking at hierarchically small scale compared to A with large anomalous dimension. 

IV. SUMMARY AND DISCUSSION 

In this paper, we studied corrections to the conformal hyperscaling relation by taking the example of SU(3) gauge 
theories with various number of fundamental fermion. From the analytical expression of the solution of the ladder SD 
equation, we identified the form of the leading correction to the hyperscaling relation. We found that the anomalous 
dimension, when identified through the hyperscaling relation neglecting these corrections (which we denoted as 7^f ), 
tends to be lower than the real value at the fixed point. 

We further studied finite-size hyperscaling relation through the ladder SD equation in a finite space-time with 
the periodic boundary condition. We found that the anomalous dimension, when identified through the finite-size 
hyperscaling relation neglecting the mass corrections as is often done in the lattice analyses, yields almost the same 
value as that in the case of the infinite space-time neglecting the mass correction, i.e., a lower value than r )* n . The 
introduction of the finite size of space-time should also break the infrared conformality, though we found that the 
correction to the hyperscaling relation due to the finite- volume effect seems to be negligible at least in the range of 
L we studied in this paper. This can be seen from the fact that the finite-size hyperscaling relation is approximately 
satisfied with 7^ which is obtained from the infinite- volume analysis. If 1/L correction were large, there must have 
been a visible violation of hyperscaling relation caused by it. The smallness of correction coming from finite-size effect 
can also be understood from the fact that a function with a form shown in Eq. (f38|) , in which only mass correction is 
taken into account, can be fitted to all the data in Fig. [5] pretty well. 

We also applied the finite-volume SD equation to the chiral-symmetry-breaking phase and found that when the 
theory is close to the critical point such that the dynamically generated mass is much smaller than the explicit breaking 
mass, the finite-size hyperscaling relation is still operative, with the mass corrections to the anomalous dimension 
being somewhat involved, however. 

From a lattice simulation point of view, there are several things we can learn from the results of the present paper. 
When the input bare mass is not small enough, and data are not precise enough to find the mass correction, finite-size 
hyperscaling plot might give fairly good aligned picture with a value of the mass anomalous dimension which is much 
smaller than the value at the IRFP. If data are precise enough, one could notice misalignment of data which is caused 
by the fact that the value of 7J™ is different for different values of the meson mass M. However, if one didn't know 
that the misalignment is fake coming from the correction term, one could draw conclusion that the theory is not 
infrared conformal, even though it actually is. As we mentioned at the end of the previous section, opposite could also 
happen, namely, even if a theory is actually in the chiral symmetry breaking phase, one could draw conclusion that it 
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is infrared conformal especially when the amount of spontaneous chiral symmetry breaking is very small and/or data 
is not precise enough. Thus, careful attitude is important when one judges whether a theory is infrared conformal or 
not by using hyperscaling analysis. However, main message of our analysis is that if one observed a certain type of the 
finite-size hyperscaling relation (with some finite mass corrections), it already hints the remnant of the IR conformal 
theory no matter it may be in the broken phase (applicable to the walking technicolor) or the conformal window: 
It implies a new situation of the 4-dimensional non-Supersymmetric gauge theories and a new phenomenological 
application. 

In this paper, we studied large Nf QCD as a concrete example for the study of hyperscaling relation, though the 
extension to different number of color and different fermion representation is straightforward. This is because, in 
the context of the SD equation with the improved ladder approximation, C2 g 2 appearing in the equation is the only 
quantity which differentiate different theories, and with the simplification of the running coupling adopted in the 
current study, this is proportional to a*/a CI . Therefore the only relevant thing is how close the value of the running 
coupling at the IRFP is to the critical coupling. 

It is also interesting to ask what is the best way of analyzing data to extract the correct picture. The SD equation 
gives us a nice playground to try to find an analysis method which works well for finding correct picture of a given 
theory since it can generate as many data as we like, and we know the "answer" , namely, whether the theory possesses 
an IRFP, and also the value of the mass anomalous dimension in that theory. We can try several different analysis 
methods with those generated data, and compare the results with the answer. By doing so, we can tell which analysis 
method produces the answer rather correctly. We tried fitting using the fit function shown in Eq. Q38[) as an example 
of such studies, and found that it works quite well extracting the true value of 7^ of the theory. Of course, it is worth 
investigating further in this direction. Various different analysis methods should be studied for the purpose of finding 
a practical method which can extract a more correct picture from lattice data. 
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FIG. 8: Values oimpL obtained from the finite- volume SD equation as a function of £ = LA(m /A) 1/(1+7) for 7 = 0,0.5, 1.0, 1.5 
and 2.0 in SU(3) gauge theory with 9 fundamental fermions. Data for LA — 12, 16, 20, 25 and 30 are plotted as different symbols. 




FIG. 9: Values of mpL obtained from the finite-volume SD equation as a function of x = LA(mo/A) 1//< - 1+7 ' for 7 = 1.0 in 
SU(3) gauge theory with 11 fundamental fermions. Data for LA = 12, 16, 20, 25 and 30 are plotted as different symbols. 



